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Suna«ry 


The MAGSAT scalar data are processed to construct a crustal magretic 
anomaly map over the continental U.S. The processing steps Include a removal 
of the reference field model, a path-by-path subtraction of a low order 
polynomial through a least-squares fit to reduce orbital offset errors, and 
a two-dimensional spectral filtering to mitigate the spectral bias 
Induced by the path-by-path orbital correction scheme. The resultant anomaly 
map shows reasonably good correlations with an acromagnetlc map derived from 
the project MAGNET. Prominent satellite magnetic anomalies are Identified in 
terms of geological provinces and age boundaries. 

To reduce the MAGSAT data to a common altitude, we attempted three types 
of functional representations of the field, namely, a spherical harmonic 
series, a finite harmonic series, and a triple half-range cosine 
series, vnille these representations properly reduce simulated satellite data 
to a common level, their results derived from the actual MAGSAT data show 
unacceptably high rms residual errors. This Implies that the difference In 
the magnetic data among different satellite paths cannot be solely accounted 
for by the variations In orbital altitude alone. Only after the temporal 
changes In the external field are properly removed, can the data be 
reduced successfully to a common level. 

Finally, we applied to the MAGSAT data a newly developed method of 
estimating the Curie depth. The Inversion method produces both the Curie 
depth topography and laterally varying magnetic suscentlblllty of the crust. 

A contoured Curie depth map thus derived shows general agreements with a 
crustal thickness map based on seismic data. 

Based on the current study, we believe that the satellite magnetometry, 
particularly with Improved methods of correcting the external field 
variation, can make further significant contributions toward the understanding 
of the earth's crustal structure. 


Analysis of the MAGSAT data over the U.S 


The HAGSAT data used for this report consist of the entire Investigator 
B data file for the satellite duration. The data area Is bounded by latitudes 
2AN-52N and longitudes 230E-295E. The magnetic field measured at the 
satellite level consists of various contributions: the Internal and external 

parts of the geomagnetic field, Its temporal variations, and the anomalous 
field caused by magnetized bodies In the crust of the earth. In order to 
study the crustal contribution of the field, we selected only those 
geomagnetically quiet time data (geomagnetic planetary Index Kp less than 2+) 
and removed the reference field of MGST (6/80) (Langel et al., 1981). 

Figure 1 shows the distribution of data as a function of satellite 
altitude. We notice that the majorlt:' (56Z) of the total 49,000 data points 
fall within the altitude range between 350-450 km. Figure 2 shows spatial 
distribution over the U.S. of about one third of the total data used for the 
study. As expected, the data are evenly distributed over the entire region. 

The wide scatter of the satellite altitude poses a rather difficult 
problem of compiling and presenting the data as a single magnetic anomaly map. 
To demonstrate this problem we show, in Figure 3, a contoured map of the raw 
data falling in a relatively small altitude window between 350-400 km without 
any corrections. By limiting the altitude range, we initially attempted to 
minimize the effect of the altitude variation. 

The orglnal data were resampled on a 1 x 1 degree grid for contouring. 

The value at each grid point was evaluated using all data within a circular 
area of a three-degree radius and imposing a weighting function of 1/R where R 
is the distance between the grid point and the data point. No reduction of 
data to a common altitude was made. 
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The orginal data thus coatoured show numerous narrow anomalies parallel 
to the satellite paths. These anomalies are believed to be caused by 
differences of the satellite altitude, static bias, and drift rate among 
paths as well as temporal variations of the external geo*nagnatlc field. 

In order to reduce these Inconsistencies among different paths, we 
removed a best-fitting low order polynominal along each path, thus minimizing 
the rms difference between the observed data and the predicted reference field 
data. The optimum order of the polynomial obviously depends on the linear 
path length and longest wavelength features desired to be retained. We have 
tested polynomials of order 0, 1, and 2. Among these, we find, somewhat 
subjectively, the first order polynominal most optimal. The resultant map, 
after removing the first order polynomials, Is shown on Figure 4. The average 
slope is found to be 3.5 x 10”^ gamma/km, with a standard deviation of 5.7 x 
10“^ gamma/km, implying that the drift rate is highly random among paths (Won 
and Son, 1982). The removal of a second order polynomial for each path did 
not change the result significantly and showed signs of over-correction 
indicated by some diminishing long wavelength anomalies. 

The removal of a best-fit polynomial along each flight path reduces 
considerably the offset problems among adjacent paths, thereby eliminating 
narrow north-south flight line anomalies. However, since the scheme does not 
take into consideration the three-dimensional nature of the magnetic field 
outside the flight paths, it can potentially create spurious long wavelength 
east-west trending anomalies. Most MAGSAT anomaly maps, including the one 
shown on Figure 4, manifest dominant east-west trending anomalies while 
suspiciously lacking lorth-soutit trending anomalies. From a mathematical 
viewpoint, once we remove a 13th order-and-degree reference field, the 
resultant data should not contain any anomalies whose spatial wavelengths are 
greater than those of the highest harmonic constituent. The fact that these 
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long wavelength anomalies still exist leads us to believe that the 
path-by-path polynomial adjustment Is the potential source of these long 
wavelength east-west anomalies. We attempted to reduce the dominant east-west 
anomalies by applying a two-dimensional spectral hlghpass filter trlth a cutoff 
wavelength of 16.5 degrees roughly corresponding to the wavelength of 13th 
order spherical harmonic term. Figure 5 shows the resultant map after the 
spectral filtering. The map Is now fairly compatible with known geological 
and geophysical features and will be the basis for geological Interpretation 
later In this report. 

Analysis of the MAGNET U.S. Aeromagnetic Data 

The MAGNET U.S. aeromagnetic survey data, obtained by the U.S. Naval 
Oceanographic Office and available from the National Geophysical and 
Solar-Terrestrial Data Center (NGSTD) of NOAA, were collected during 1976 and 
1977 at altitudes above ground level of 500-700 m In non-mountalnous terrain 
and 900 -1,000 m in mountainous terrain. The MAGNET file contains some 
650,000 measurements at approximately 100 m Intervals along the flight path. 
Spacing between the N-S flight paths is about 1.1 degree in longitude. Since 
the data density along the path Is much higher than the line spacing, we 
averaged over every 100 data points to reduce the total data amount to about 
6,500 points. This gives approximately a 0.1 degree sampling interval 
along the path. 

Initially, we removed a I3th order-and-degree reference field of 
GSFC(9/80-2) with secular variation. The residual field thus obtained still 
showed considerably longer wavelengths than those of the 13th order spherical 
harmonics. Although the reason for this discrepancy Is not clear, It Is 
perhaps due to Insufficient removal of the non-crustal field. Other reference 

field models produced varying, yet similar results. 
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In order to reduce further the remaining long wavelength anomallea from 
the data, we removed a beat-fitting third order polynomial surface from the 
MAGNET residual data. The resultant data were contoured from 1x1 degree 
gridded data with the same sampling %rindow and weighting function used for the 
HAGSAT data. The map (Figure 6) thus produced shows a qualitative resemblence 
to the MAGSAT anomaly map of Figure 3. 

Crustal Geological Interpretation of MAGSAT Data 

The MAGSAT map (Figure 7) shows mostly broad-scale features that reflect 
fundamental crustal structures that persist to depths at which the Curie point 
is reached. Although some upper crustal structures can be delineated on the 
map, the finer details that can be correlated on standard aeromagnetic maps 
(Zietz, 1980) are not present. The anomalies on the MAGSAT map must be 
related to structures that exist In the lower crust and possibly even the 
upper mantle In some areas. 

Northeast Trending Ano mal ies 

Perhaps the most striking feature of the MAGSAT crustal anomaly map is 
the northeast trend of the major magnetic anomalies (Figure 7). Northeastern 
trending structures have been described across the entire United States and 
represent one of the principal crustal fabrics of the crayon. Some of these 
better documented structural features are: 

(1) Colorado lineament - a belt of Precambrlan faults (Warner, 1978) 

that cuts across the Colorado Plateau from south of the Giand Canyon 
and Rocky Mountains into southeastern Wyoming (Figure 7). It Is 
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approximately allftoed with the boundary of two distinct Precambrlan 
terrains In Minnesota (Morey and SimSt 1976). 

(2) Proterozoic chronologic province boundary that separates 1690-1780 
Ma old rocks on the north from 1610-1680 Ma old rocks on the south 
(Silver, 1968; Silver, et al., 1977). This boundary extends 

from southwestern Arizona through northern New Mexico Into 
southeastern Colorado. 

(3) Shear zones In the western United States (Chapin, et al., 1978) and 
In the Great Lakes region (LaBerge, 1972; Sims, 1976). 

(4) Midcontinent geophysical anomaly that extends from Lake Superior to 
Nebraska and Is considered to be a major crustal rift. It contains 
Keweenawan basalts and Is flanked by elongated basins filled with 
red sandstones (King and Zletz, 1971). 

(5) New York“Alabama lineament which Is defined by a series of linear 
magnetic anomalies and separates provinces with different trending 
gravity anomalies vKing and Zletz, 1978). Zletz (1980) has pointed 
out that the lineament seems to delineate the northwestern boundary 
of a selsmlcally active block In the eastern United States. 

(6) Selsmlcally active belt that extends from the St. Lawrence River 
Valley through New Madrid, Missouri (King and Zletz, 1978). 

It Is clear that northeast-trending structural features form some of the 
largest and most pervasive disruptions in the crust across much of the United 
States. They all pass through regions that overlie Precambrlan basement rocks 
at depth and are well Inboar'd of tnose areas thought to represent suspect 
terrains. The parallelism of the northeast-trending crustal structures and 
the major NAGSAT linear anomalies suggests the two are directly related. 
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Prtcaabclan Terrains 


Prccambrisn terrains also show soae correlation with the MAGSAT data. 

One of the aost striking correlations Is the apparent coincidence of the 
southern boundary of Archean gnelss-islgBatlte baseoent (Van Schaus and 
Bickford, 1981) and the aagnetlc low that sweeps across the aldcontlnent froa 
Wisconsin to Colorado and then northwestward Into Idaho. This boundary Is 
thought to aark the edge of the Archean craton and, therefore, aay represent 
the site of ancient plate tectonic processes. 

Another, less convincing, correlation is the approximate coincidences of 
the southern boundary of the early Proterozoic (1690'~1780 Ha old) belt of 
aetasedlmentary, metavolcanlc, and plutonlc rocks that are largely felslc (Van 
Schraus and Bickford, 1981) and the northern limit of the belt of magnetic 
highs that extends from Texas to Michigan. There Is no agreement as to the 
tectonic character of this belt of rocks; both plate margin and Intracratonlc 
settings have been proposed to explain Its origin. Van Schmus and Bickford 
(1981) argue rather convincingly that the available data Indicate the be^t 
represents a continental margin. 

Geological Provinces 

On a smaller scale, correlations can also be made with certain geological 
provinces. Some of the more obvious are: 

(1) The Northern Rocky Mountains are magnetically high. 

(2) The Middle Rocky Mountains are generally low magnetically except for 
a weak magnetic high that parallels the Front Range. 

(3) The Great Basin of the Basin and Range Provinces Is characterized by 
a magnetic low that may be the result of decrease In depth to the 
Curie point In that region. 
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Rdductioa of Satolllt* 0«t« to a Gowaon AltlCud* Surfaco 


In ordar to ranova tha affect of altitude variation and to adequately 
praaarve broad-acala ugnetlc anoaaliaa, m next attaapt to reduce the data on 
a conaon altitude. Thera are two principal aethoda: 1) equivalent aourca 
representation and il) functional representation. 

Mayhew (1979) uACd the equivalent source techniques where the paraaetars 
for regularly-spaced aagnetlc dipoles are deteralned in the least-squares 
sense so chat the aagnetlc field calculated froa such a source distribution 
approximates the observed data. Bhattacharyya (1977) proposed a technique 
specifically designed for the efficient reduction of satellite data to a 
conoton altitude and inclination. From an equivalent magnetic dipole 
distribution of magnetization at surface, he calculated the magnetic field at 
any point above the surface utilizing the relationship between the observed 
values and the spatial gradient of magnetization. 

As for the functional representation approach, we have considered a) a 
spherical hariK>nlc series, b) a finite harmonic series in Cartesian 
coordinates, and c) a triple half-range cosine series In Cartesian coordinates 
(Son and Won, 1982). 

First we subjected the entire data over the U.S. to a least-squares fit 
to the surface spherical harmonic series up to 15th order-and-degree. The 
resultant coefficients are then used to recompute the field at a common 
..iCltude of 350 ka (Figure 8). The result Is rather disappointing: while the 

re> iltant map retains the broad features. It Is completely devoid of small 
features of geological Interest. The overall rms error amounts to about 5 
gammas. Due to reasons to be discussed later, a further Increase in the 
number of terms does not Improve the resolution; nor would the reduction of 
area size resolve small wavelength features. 
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N«xc, «• consider the finlts hsraonic ssrlss rsprsssntstlon* Henderson 
end Cordell (1971) used this rspresentstion coabincd with the lesst-equaree 
technique for three dlaenslonsl data. Regan (1979) applied It to reducing 
satellite eagnetlc data using Cartesian coordinates. We applied this aethod 
to the MACSAT data over a portion of the eastern U.S. 

Figure 9 shows a contour nap of raw data In the region uncorrected for 
orbital offset errors. Figure 10 shows an anomaly map recomputed at 400 ka 
altitude by the finite harmonic series representation using 10 terms along 
each direction (100 coefficients). The overall rus error Is still 2.9 gammas 
unacceptably high. 

Finally, we consider a triple half'*range cosine series In the 
Cartesian coordinates to represent the crustal magnetic field T(x,y,t): 

«o »v> lo * y * 

T(x,y,r) "I I I cosr m - cosx n - cost 1 - (1) 

m*0 n"0 1*0 Ly 

where L^, Ly, and are the dimensions of an arbitrary rectangular volume 
occupying the data space. The coefficient can be found by the method of 

least-squares. Once these coefficients are determined for the data they may be 
used to recompute the field at any point inside the volume. The number of 
terms to be uaeJ in the series approximation mo, Oq, 1q depends on the size of 
the region to be analyzed and the required spectral resolution of the crustal 
field. The method has been tested successfully against sets of simulated 
satellite data due to simple anomalies. Figure 11 shows the resultant anomaly 
map through the triple half-range cosine series representation using 6 terms 
along each direction (216 coefficients) in the eastern U.S. region. The 
overall rms error Is 2.4 gammas. 

From all three methods and results, we do not note any significant 
Improve iients on the anomaly resolution through the reduction to a common 
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altitude* Let ue 'roatider lome beeic reeeone for this apperent failure ia the 
followiog. 

All three repreaentatioaa are the aolutioaa of the Laplace equation in 

the respective coordinates and, therefore, physically valid* It should be 

noted, however, that the raw data subjected to these representations were 

inevitably not corrected for orbital offsets* Thus, if these orbital offsets 

are truly caused by the altitude variation alone, the oethods should have 

« 

properly reduced the data to a common altitude to the spatial resolution of 
the highest order base function used for the field representation. However, 
it is obvious that these results do not show the spatial resolution prescribed 
by the series representation. 

This fact, along with high overall rms errors, testifies that the 
altitude variation alone cannot completely account for the orbital offsets. 

It appears that other causes such as the temporal variations of the external 
field dominate the orbital offsets. Only after these causes are properly 
accounted for and removed, the representation method can successfully reduce 
the data to a common altitude. Presently, ve do not seem to have solution for 
this difficult problem. 

Interpretation of MAGSAT Data for Crustal Magnetization 

The basic objective in relating the magnetic field data %rith the crustal 
structure of the t.'^rth is to compute the lower depth limit of magnetized 
masses in the earth's crust. Rocks lose their magnetism at the Curie 
temperature at which ferrlmagnetic rocks become paramagnetic, and their 
ability to produce detectable magneclzation disappears. Thus, the deepest 
level in the crust containing materials with discernible magnetization is 
generally Interpreted as the depth to the Curie point isotherm. 
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The Curie point le about S80*C for aegnetite. With epproprlete titenlun 
•ubetltutiona, Buddlogtoo end Lindeley (1964) calculated an /average Curie 
point ranging between S20*C and 560*C for rocka in the deep It ia 

generally believed that the aaount of geotheraal heatflow should correlate 
with the Curie depth and thus, in turn, to the crustal aagnetic field. 

Our aain goal is, therefore, to determine the bottom shape of the 
magnetised crust from a magnetic anomaly map. Since the magnetic anomalies 
attributable to the bottom geometry are usually quite smaller and have much 
longer wavelengths than those produced by shallow geological variations, the 
problem is comparable to searching for a needle in a haystack. Early studies 
include those by Vacquler and Affleck (1941) and Bhattacharyya and Morley 
(1965). In both cases, each isolated anomaly was filtered and separately 
interpreted by the empirical graphic method using a vertical-sided prism. 

A more sophisticated method was proposed by Bhattacharyya and Leu (1975a 
and 1975b). Their method requires an extensive initial filtering of the 
magnetic data in both regional and short wavelength domains. The filtered 
data Is subsequently divided into a large number of blocks. For each block, a 
two-dimensional spectrum and its moments are computed and compared with a 
model of an isolated vertical-sided prism within a block in order to locate 
the corners of the body. The total amount of computation is tremendous since 
the method requires a two-dlnensional Fourier transform for each block. 
Applying the method to the Yellowstone National Park area, they produced the 
Curie isotherm map well correlated with the known geothermal area. 

Employing a similar technique, Shuey et al. (1977) concluded that it is 
essentially impossible to determine Curie depth with any resolution at all by 
fitting a vertical prism to a single anomaly. The Curie depths they derived 
could be changed by as much as 10 kn without violating the observed data. 
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This conclusion Is seeningly in conflict vlth those of Bhstcschsryys end Leu 
(1975b). 

All aethods reviewed here ere coaaonly besed on the essuaptlon that there 
exists an Isolated aagnetlc source for each anomaly. Each Individual anomaly 
is assumed to be caused by a single vertlcal^slded prism (Bhattacharyya and 
Leu» 1975a, and 1975b) or a truncated vertical cone (Shuey et al., 1977). 

Such Isolated models are apt to generate 'Spurious anoinalles due to their 
unrealistically vell~dsfined corners and vertical surfaces. These spurious 
anomalies can induce significant errors in either direct-modelling or spectrum 
c«lcui.ation. 

Rock formations causing long wavelength laagnetic anomalies at a depth 
close to the Curie point are more likely to have a continuous lateral 
dlstritution rather than isolated blocks of well-defined geometrical bodies. 

A realistic model at this depth should manifest a continuous lateral 
distribution of magnetic materials havinp, variable thicknesses and 
susceptibilities. 

Fluctuations in long wavelength magnetic anomalies can be attributed to 
lateral variations either in magnetization strength or in Curie depth. These 
double uncertainties make the task of simultaneously determining both the 
magnetizations and the Curie depths very difficult, if not impossible. 

Similar uncertainties apply to many geophysical modelling theories, e.g>, a 
thin magnetic dike for which the anomaly Is the same, as long as the product 
of thickness and susceptibility remains the same. However, it can be shown 
that the statement is no longer true if the dike has a considerable thickness, 
for which case both the thickness and the susceptibility can be independ ntly 
determined from observed data (Won, 1981). The r*eient approach is based on 
the classics? Gauss method for solving non-linear equations (Corbato, 1965; 
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Johnson^ 1969; Won, 1981) coupled vlth Merquerdt'e inversion Mthod 
(Merquardt, 1963) to derive continuous crustel thickness end susceptibility 
profiles from regional magnetic date (Von, 1982). 

Figure 12 shows the model which is used for inverting msgnetic date. The 
model consists of laminated thick vertical prisma having flat top surfaces and 
linearly-connected inclined bott<» surfaces. The magnetic susceptibility 
below the lower boundary is assumed to be zero so that the bottom geometry 
represents the Curie isotherm topography. Althcagh data will be confined 
within the laminated block region, two seBi-lnf Inlte slabs are added on either 
side in order to reduce the edge effects of the first and last blocks. The 
unknown parameters to be determined are the depth (h*s) at each nodal point 
and the magnetic susceptibility (k's) of each prismatic body. 

The model Is two-dimensional with an arbitrary strike angle with respect 
to the magnetic north. Since the method uses total field mag'>etlc data, there 
is no need for reducing the data to the polar anomalies. The magnetic anomaly 
generated by a single vertical block having a flat top and an Inclined bottom 
can be derived by analytically combining two Inclined dikes. By summing up 
these Individual blocks, we can compute the total field for any given set of 
blocks having variable depths and susceptibilities (Won, 1982). 

Techniques for determining unknown parameters of a nonlinear function 
Involve Iteratively correcting currently assumed parameters by differential 
amounts, thereby minimizing the rms error between the theoretical prediction 
and the observed data. Two predominant techniques of determining the 
.orrectlon amounts are Gauss' method and the gradient, or the steepest-descent 
T'ethod. Marquardt's method combines these methods by controlling the amounts 
of ditterentlal correction to Insure both convergence and speed. 
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Using th« geosetricsl sodel sod the iavsrsion tschaique thus described, 
ve snslyted s MAGSAT crustal snoosly profile scross U.S. coutineat st leittude 
36*N. Figure 13 shows the dste profile end the computed profile slong with 
the bottosi topography of the magnetised crust and the lateral susceptibility 
variation. The ms deviation of the model is only 0.19 gammas. 

A total of twelve equally-spaced east-west MAGSAT profiles were subjected 
to this inversion technique and n Curie depth map is constructed (Figure 14). 
As expected, the magnetized crust is thin in the oceanic regions (less than 30 
km) with relatively high magnetic susceptibilities (see Figure 13), while it 
is very deep over the geologically old Apalachian region. Comparison of this 
map with the crustal thickness map (Figure IS) compiled by Allenby (1980) from 
seismic data shows qualitative agreements. We believe that the inversion 
technique can be improved even further by adopting a three-dimensional crustal 
magnetization model. 
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1. MikGSAT scalar data population as a function of satellite altitude 
over the continental U«S. 

2. Footprints of one-third of MAGSAT data In the continental U.S. 

3. MAGSAT raw residual data within an altitude window of 350-400 ka 
contoured at a two gaoma Interval. 

4. MAGSAT anomaly s^sp after removing a best-fit straight line from 
each path* 

5. MAGSAT anomaly map after a two-dimensional highpass spectral 
filtering for wavelengths shorter than 16.5 degrees. 

6. MAGNET U.S. aeromagnetlc map referenced to GSFC (9/80-2) field 
model contoured at a 50 gamma Interval (after Won and Son, 1982). 

7. Geological correlation of the filtered MAGSAT crustal anomaly map. 

8. Reduction of MAGSAT data to 350 km altitude by a spherical harmonic 
analysis contoured at a two gamma Interval. 

9. MAGSAT raw data over a portion of the eastern U.S contoured at a 
two gamma Interval. 

10. An anomaly map of the eastern U.S. reduced to 400 ka altitude 
through the finite harmonic series representation. 

11. An anomaly map o^ the eastern U.S. reduced to 400 ka altitude 
through the half-range cosine series representation. 

12. Model geometry used for inversion of magnetic data to estimate the 
Curie depth and lateral susceptibility variation. 

13. Profiles of the bottom of the magnetized crust and susceptibility 
based on a MAGSAT crustal anomaly profile at 36*N across the 
continental U.S. 

14. Thickness of the magnetized crust based in the continental U.S. 
on the MAGSAT data. 


Figure 15. Crustal thickness based on seismic studies (from Allenby, 1980). 
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Figure I. MAGSAT scalar data population as a function of satellite altitude 
over the continental U.S. 
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MAGSAT DATA POINTS ( I \ 1 1173 - 5126180 , Kp < 2 ) 
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Figure 2. Footprints of one-third of MAGSAT data In the continental U.S. 
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Figure 3> MAGSAT raw residual data within an altitude window of 350-400 ka 
contoured at a two gana Interval* 



ORIGINAL PAGE IS 
OF POOR QUALITY 



Figure 4. MAGSAT anonaly nap after removing a 
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Figure 5 


MAGSAT anonaly map after a two-dlaens tonal hlghpass spectral 
filtering for wavelengths shorter than 16.5 degrees. 
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Figure 9 
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MAGSAT raw data over a portion of the eastern U.S contoured at a 
two gamma interval* 
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Figure 10. Aa anomaly map of the eastern U.S. reduced to 400 km altitude 
through the finite harmonic series representation. 
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Fl){ure 12. Model geometry used for Inversion of magnetic data to estimate the 
Curie depth and lateral susceptibility variation. 
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Figure 13. Profiles of the bottoa of the aagnetlzed crust and susceptibility 
based on a HACSAT crustal anonuily profile at 36*N across the 
cuuitneaLal U.S. 
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Figure 15. Crustal thickness based 


on selsaic studies (from Allenby, 1980) 





